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We, Gabor Forgacs, Karoly Jakab, Adrian Neagu, and Vladimir Mironov, declare 
as follows: 

1 . We are the co-inventors of the subject matter claimed in the above- 
referenced patent application, U.S. Patent Application Serial No. 10/590,466 (the '466 
application). 

2. We conceived of and reduced to practice the subject matter of the 
invention as claimed in pending independent claim 52 prior to October 4, 2003. Exhibits 
A and B attached hereto are submitted as evidence of our conception and reduction to 
practice prior to October 4, 2003. 

3. Exhibit A is a true and correct copy of an article by Jakab et al., which was 
published in the journal Proceedings of the National Academy of Sciences (PNAS) on 
March 2, 2004. The face of the article indicates that a manuscript for the article was 
received for review by the journal on July 20, 2003. We are all co-authors of this article, 
and conceived of and reduced to practice the subject matter described therein. 

4. Exhibit B is a true and correct copy of the manuscript for the Jakab et al. 
article originally submitted to the journal PNAS. It will be appreciated that the 
manuscript as originally submitted to the journal PNAS is substantially similar to the 
published Jakab et al. article, and that both documents describe the subject matter of 
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independent claim 52 of the '466 application.^ In particular, the manuscript describes: a 
three-dimensional layered structure comprising at least one layer of a matrix; and a 
plurality of cell aggregates, each cell aggregate comprising a plurality of living cells; 
wherein the cell aggregates are embedded in the at least one layer of matrix in a non- 
random predetermined pattern, the cell aggregates having predetermined positions in 
the pattern (see, e.g., p. 2, lines 16-17 & 21-23; p. 4, line 21-p. 6, line 18; p. 14, lines 
7-13, p. 25, lines 1-11; and Figure 3). 

5. Based on the foregoing, we believe that we conceived of and reduced to 
practice the claimed invention prior to October 4, 2003. 

6. All work referred to herein was carried out in the United States. 

7. We further declare that all statements made herein of our own knowledge 
are true and that all statements made on information and belief are believed to be true; 
and further, that these statements were made with the knowledge that willful false 
statements and the like so made are punishable by fine or imprisonment, or both, under 
18 U.S.C. § 1001, and that such willful false statements may jeopardize the validity of 
the application or any patent issuing thereon. 



^ Figure 1 of the manuscript as originally filed could not be located, but was substantially similar to Figure 
1 of tlie Jal^ab et al. article as publislied, the primary difference being that Figure 1 as published included 
an additional collagen concentration (1.2 mg/ml). This is clear from a comparison of the figure legends 
for Figure 1 in the manuscript (page 25) and the published article (page 2866). 
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Communicated by Joe! L. Lebowitz, Rutgers, The State University of New Jersey, Piscataway, NJ, January 8, 2004 (received for review Ji.iiy 20, 2003) 



Seif-assembiy is a fiindamentai process that drives structural or- 
ganlzatiors im both jnariimate and living systerns. it is in the course 
of seif-assembiy of ceSis and tissues m easiy deveioptrtent that the 
organism and its parts everstualiy acquire their finai shape. Even 
though deveioprrsenta! pattemirsg through se!f-assembiy is under 
strict genetic control it is cisar that uStlmately it is physical mech- 
anisms that bring about the complex structures. Here we show, 
both experimentally and by using computer simulations, how 
tissue liquidity can be used to buiid tissue constructs of prescribed 
geometiy in vitio. Spherical aggregates containing many tiioii- 
sands of ceils, which form because of tissue Siquidity, were im- 
planted contiguously into biocompatible hydrogsis in circular ge- 
ometry. Depending on the properties of the gei, upon incubation, 
the aggregates either fused into a toroidal 3D structure or their 
constituent cells dispersed into the surrounding matrix. The model 
sirnuiations, which reproduced the experimentally observed 
shapes, indicate that the control parameter of structure evolution 
is the aggregate-gel interfacial tension. The model-based analysis 
also revealed that the observed toroidal structure represents a 
metastable state of the cellular system, whose lifetime depends on 
the magnitude of celi-ceil and cell-matrix interactions. Thus, these 
constructs can be made long-lived. We suggest that spherical 
aggregates corrsposed of orgarv specific cells may be used as 
"bio-ink" in the evolving technology of organ printing. 

Self-assembly is the fundamental process, which generates 
strtictura! organization across scales (1). Histogenesis and 

org;iiu:'gciK:si:s examples of self-assetribly processes, in 
wliicli, lliiough ccli-cell and ceii-extracelkuai matrix interac- 
tions, the developiiig oigaiiisrfi and its parts gradually acquire 
their final shape. In the preset! t worli we use both experimental 
and computational approaches to demonstrate how the self- 
organizing properties of ceils and tissues, the basis for morpho- 
genesis, can be exi)i()ited to build 3D biological structures of 
prescribed geojneti-y. 

Tissue engineering (2-7) offers the opportunity to study 
self-assembly processes diiring histo- and organogenesis in vitro, 
under controlled conditions. Tissue engineet ing aims not only to 
create desirable organs, but also to better understand the 
fundamental mechanisms and principles of biological organiza- 
tion in general and morphogenesis in particular. Clas.sical tissue 
engineering is based on seeding cells into biodegradable polymer 
scaffolds or gels, cuituring and expanding thern in bioreactors for 
several weeks, and finally implanting the resulting tissue into the 
recipient organism, where the maturation of the new organ takes 
place. 

It has recently been suggested to use cell aggregates, instead 
of itidividuai cells, as building blocks in tis.sue engineering (8, 9). 

Cell aggregates have traditionally been used as a powerful tool 
to understand the principles of cell-ceil (10) and cell-matrix 
adhesion (11), as well as ceil sorting (12). In addition, rapid 
prototy[>ri!g technology lias successfully been applied for com- 
puter-aided deposition of ceils in gels to create 3D tissue 
constructs (13, 14). We suggest that cell aggregates may be used 
as drops of "bio-ink," which, upon implantation or "printing" 



into a scaffold ("bio-paper"), have the ability to fuse into 3D 
orgai! structures (15-17). 

The ability of cell aggregates to fuse is based on the concept 
of tissue fluidity (18, 19), according to which embryonic tissues 
in many respects can be considered as liquids, fn particular, in 
suspen.sion or on nonadhesive .surfaces, varioiis mtiiticelhilar 
aggregates round tip into spherical shape similarly to liquid 
ci'-opJets (12). We hypothesize that closely placed aggregates in 
.:pproprialely chosen 3D gels can fuse to form tissue constructs 
of desired geometry. 

To demonstrate the feasibility of such a proposition, we used 
aggregates of genetically transformed cells with controlled ad- 
hesive properties, and arranged them to form a ring in gels of 
different chemical and mechanical properties. Our results deir;- 
onstrate that contiguous aggregates under appropriate condi- 
tions, defined by the composition of the embedding gel, indeed 
can fuse into structures of specified morphology. We have also 
constructed and experittientally validated a mathematical model 
of cell aggregate fusion. 

Tciken together, these j'esuits suggest that cell aggregates can 
succes.sfuHy be used as building blocks it5 tissue engineering 
technologies. They iilso provide support for the concept 
self-assembling bio-ink, and thus justification for lire use of cell 
aggregates in the evolving organ printing tools. 

filateriak aiid ^Vlethods 

Ceil Aggregate Preparation. Chinese Hamster Ovary (CHO) cells, 
transfected with .iV-cadherin (couttesy of A. Betshadsky, Weiz- 
mann Institute, Rehovot, Israel), were infected with histone 
binding H2B-YFP retrovirus (kindly provided by R. D. I.ans- 
ford, Beckinari Institute at California Institute of Technology). 
Confltient cell cuhxires (3-4 X 10* ceils per 7S-cm^- TC dish) 
grown ifi DMEM (GIBCO/BRL) suppiemertted with 10% FBS 
(U.S. Bio-Technologies, Pottstown, PA), 10 fig/mi penicillin, 
streptomycin, gentamicin, and kanatnycin, 400 fxg/ml gene1;i- 
cin), were washed twice with Hanks' balanced salt soiulion 
(HBSS) containing 2 mM CaCla, then treated for 10 min with 
trypsin 0.1% (diluted from 2.5% stock, GIBCO/BRL). De- 
pleted cells were centrifuged at 2,500 rpm for 4 min (Fisher 
Centrific model 225). The resulting pellet was transferred into 
capillary micropipettes of 500-jOi.m diameter and incubated at 
37''C with 5% CO2 for 10 min. The firm cylinders of cells 
removed from the pipettes were cut into fragments (500-fj.m 
height), then incubated in IQ-mi tissue culture flasks (Bellco 
Glass, Vineland, N.T) with 3 ml of DMEM on a gyratory shaker 
at 120 rpm wjth'5% CO? at 37*C for 24-36 h. This procedure 
reproducifaly provides spherical aggregates of simrlai srze 
('■■«500-pitn diameter). 



Abbreviations: Cf-iO. Chinese hamster ovary; MCS, Monte Cario stc-p; DAH, differential 
adhesion hypothesis. 

^To whom correspondence should be addressed, t-mail: forgacsgS'.'nissouri.edL;. 
© 2004 by The National Academy of Sciences of the USA 
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Ceii Aggregate-Gel Stsructures. NeuroGel (a biocompatible poroiis 
poly[A'-(2-iiydroxypropyl)methacrylamide] hydrogel) disks of 
10-mm diameter sine! 2-Enfn thickness, containing Arg-Giy-Asp 
(RGD) fragments (kindly provided by Stephane Woerly, Or- 
ganogel Canada, Quebec) were washed three times witli DMEM 
to eliminate ihe storage medium. This gel has been shown to 
provide favorable condiiions for spinal cord repair (20, 21). A 
0..5-n!rt! wide, 0.5-ssin). deep circular groove was cut into a disk, 
then filled with 10 aggregates, placed contiguously to form a 
closed circle. The groove was refilled with the gel to completely 
embed the aggregates. This strticture was incubated at 37°C, 5% 
CO2 for 72 h in a tissue culture dish containing 10 m 1 of DMEM, 
washed with PBS, and finally embedded in Tissue-Tek OCT 
Com.pound (Electron Microscopy Sciences, Fort Washington, 
PA). The structure was slowly cooled {1°C per min) to — 20°C in 
a Nalgene freezmg container (Naigene Labware, Rochester, 
NY ). To visualize aggregate fusion at the end of the experiment, 
cryosectioning was performed with a Reichert 2800N Frigocut 
cryotome (Reichert Jung, Arnsberg. Germany), yst'lding 10- to 
16-fj.in-thin slices mounted on tnicroscope slides. Slices were 
examitjed on an Olympus IX-70 inverted microscope with flu- 
orescent attachment at X4 magnification. 

To tune the strength of cell-gel interaction, further fusion 
experiments were conducted in rat-tail collagen type I (Sigma). 
Collagen was dissolved in 1 M acetic acid, then treated with 
Ham's F12 medium with sodium bicarbonate. At room temper- 
ature, this mixture gels in a few minutes depending on concen- 
tration. The gel-aggregate structure was achieved by creating a 
ring of 10 aggregates on ilie top of a previously (alniost) 
solidified collagen layer, then covering with liquid coilageii that 
eirihcdJoJ the aggregates after gelation. These samples were 
incuixiied under the same condiiions as described above. Work- 
ing with i.O, 1.2, and 1.7 mg/ml collagen, the samples were 
II :iii:sparciit; thus, it wns ptissiblo to follow pattern evoltitiot) in 
time by piuihc coiiliasi and fluorescent microscopy. 

Cell survivability was checked with Trypan Blue (GIBCO/ 
BRT) at the end of each fusion ex|)eriment. A minimal number 
of uniformly distributed dead cells were found. 

iViodeiiiig Structurs f oimation by IVieans of Ceil Aggregate Fusion. To 

investigate sliapc cliariges c-f ihe evolving pattern, we constructed 
a .si tuple 3D inudel, it! which tlic sites of a cubic lattice are 
occupied either by point-like cells or gel volume elements. The 
total interaction energy, E of the system is written as 

where /■ and r' label lattice .site.s, and {r, /■') .signifies sumtriation 
over neighboring sites, each pair counted once. First, secoitd, and 
third nearest neighbor are included, and we assume that a cell 
interacts with the same strength with all of the 26 cells it comes 
into contact with. To specify occupancy, we assign a spin value, 
<r, to each lattice .site with values 0 for a "gel particle" and 1 for 
a cell. The interaction energy of two neighbois, /(a,., o-, ), ntay 
take either of the values i(0, 0) - -Bgg, /(I, 1) - -- s^c, or J (6, 
1) = /(I, 0) = "Ecg- Here the positive parameters egg, ecc, and 
Beg accoiiiit foe contact interaction strengths for cell-cell, gel- 
gei, and cell-gel pairs, respectively. More specifically, these are 
mechanical works needed to disrupt the corresponding bonds. 
(Note that eco and egg are works of cohesion, whereas s^g is work 
of ailhesion per bond; ref. 22.) The strength of cell-cell inter- 
action may be deterttiined experimentally either directly (23) or 
by iTicasuring the tissue surface tension (12, 24). The cell-gel 
interaction is tunable via the concentration of RGD groups in 
the gel (6, 25) or by the concentration of collagen. The gel-gel 
"bond energy" Ls an effective measure of gel filament density. 



interaction, and stiffness. It is determined by the specific cheiri- 
istry of the gel. 

The energy in Eq. 1 may be rewritten by separating interfaciai 
and bulk terms in the sum. As a result we obtain 

ycgScg + const. [2] 

Here, Bng is the total number of cell-gel bonds, and Vcg =^ (scc 
+ «!gg)/2 - Ccg is [)roportional to the cell-gel ititerfacial tetssion 
(2.2, 26, 27). The remaining terms in E do not change as the 
cellular pattertt evolves. Otir mode! is inspired by earlier efforts 
aiming at computer simulations of cell sorting (26-28), the 
morphogenetic phenomenon in which one of two, initially 
randomly interttiixed cell populations sorts out and becomes 
surrounded by the other. 

The evolution of the system is followed by using Monte Carlo 
simulations (29), relying on a randotn number generator of 
L'Ecuyer with Baym-Durham shuffle (30). The program iden- 
tifies the cells on ihe aggregate --gel interface, picks otte of tlietn 
randomly, and exchanges it with an adjacent gel particle chosen 
by chance. The corresponding change in adhesive energy, AE, is 
calculated, and the new configuration is accepted with a prob- 
ability P = 1 if < 0 01 P = exp(-j3A£:) if AE > 0. ^ = l/Ej 
is the itiverse of the average biological fluctuation energy Ej, 
analogous to the thermal fluctuation energy {31), kuT (where k-a 
is Boltzmann's constant and T is the absolute temperattire). In 
statistical mechanics, this enetgy characterizes thermal agitation 
in a system of atoms or molecules. In the case of cells, it is a 
measure of the spontaneous, cytoskeleton-driven motion of 
cells, able to break adhesive bonds between neighbors via 
membrane raffling (32), or more generally, via membrane 
protrusive activity (e.g., filopodial extensions). By definition, a 
Monte Carlo step (MCS) or ''unit oi' time," is completed when 
each cell in contact with ihe gel has iKcn given the chance to 
move once. During each MCS, the into; lacial sites are selected 
iii ratidom order. The gel boundary is ttealcd as a fixed physical 
limit of the system, and cells are constrained to move within the 
gel. 

Results 

Structure Formation: Sirsiulations, The resiilt stated in Eq. 2 indi- 
cates that the evolution of the cellular pattern is governed by a 
single parameter, ycJEj, which, for cells with specific adhesion 
apparatus, is controlled by the chemistry of the gel. The theo- 
retical analysis shows that, once gel properties are appropriately 
tutted, efficient fusion of adjacent aggregates takes place. This is 
illustrated in the simulation shown in Fig. 1. For small y^JEj (- 
0.25 in Fig. 1 K and L), cells can spread in the btilk of the gel (i.e., 
permissive gel) aiici the pattetn evolves toward its lowest energy 
state, being a sphere. (Because the interfaciai energy is small in 
comparison to the fluctuation energy, the .spheroidal structxtre in 
Fig. IL is rather disperse.) Under optimal cell-gel interface 
properties, expressed in our model by a certain range in ycJEi, 
fusion of aggregates results in a 3D toroidal structure. An 
example of this kind is depicted in Fig. lA andfi {y^JEi = 0.9, 
nonperraissive gel). 

Energetics Ourtng Structure Formation. Some shapes cijrrespotsd f:o 
local minima of the interaction energy. These represent meta- 
stable configurations. They are identified from plateaus in the 
plot of the total interaction energy vs. MCS, and are important 

for tissue engineering, for they can be m.ade long-lived. This is 
illustrated in more detail in the simtilation shovvn in Fig, 2, where 
the initial state progresses toward a metastable toroidal config- 
uration, whose energy is essentially unchatsged in the eiitire 
interval between 10"* and 6 X 10* MCS (Fig. lA). Eventually the 
toroid becomes unstable, and at <^ \aP MCS it ruptures (Fig. IB). 
Subsequent massi%'e rearrangements lead to a pronounced en- 
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Fig, 1, Initial (Upper) and final (Lower) ceii aggregate configurations in tlie simulations and in experiments using various biocompatible gels. Panels A and 6 
correspond to simulations with -/cg/f r - 0.9, and K and L correspond to simulations with ycg/Ej - 0.25. The 10 aggregates, each containing 925 cells, are one 
cell diameter from each other in the starting configurations. The final configurations are reached after 25,000 and 50,000 MCS, respectively. C-J correspond to 
CHO cell aggregates embedded in a neurogel with RGD fragments (Cand D) and collagen gels of concentration 1 .0 (f and F), 1 .2 (G and H), and 1 .7 (/and J) mg/ml. 
The nuclei of the cells are fiuorescently labeled (see Materials and Methods), and the images of the cellular patterns were acquired with a .x4 objective. The 
average diameter of the aggregates is 500 /xm. 



ergy ciecfease while the systeEn evolves inl.o ihree rounded 
aggregates. The.se remain stable lor a. long iirne liecause large 
spatial separaticMisi hinder their fusion inlo a single spheroid. 
Similar simulations showed that metastabUity depends on both 
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Fig. 2. (A) Metastabie configuration at Tcg/EV =1,1, Fusion of 16 contiguous 
aggregates (1 23 cells each) takes place within the first 5,000 MCS. The plateau of 
the energy vs. MCS corresponds to a torus. (6) TJie evolution of the energy during 
5 X 10^ MCS and representative shapes. The final, rounded aggregates remained 
stable, and the energy was practically constant up to 10- MCS (not shown). 



system size and interactiop. strengths. Because the evolution of 
the cellular pattern is driven exclusively by energy miniraization, 
and "time" is measured la MCS, in its present form, the model 
casinot provide information on the true dynamical behavior of 
the system. 

Once the structure reaches the metastabie state, it can be 
stabilized by dissolving the supporting gel. In the simulations, 
this corresponds to increasing the value of jcg/Ej. Indeed, if in 
the simulation shown in Fig. 2A this quantity is changed to 
Ycg/'^T — 2 anywliere in the plateau region, the energy remains 
constant as long as the simulation is run (results not shown). 

Experimenta! Realization of 3D Structure. To study the feasibility of 
engineering 3D tissue constructs of prescribed geometry, we 
have "manually printed" (i.e., embedded) aggregates of living 
cells into biocompatible gels. As indicated by our model, the 
ability of aggregates to fuse depends on the mutual properties of 
the cells and gel. The results of our experiments in Fig. 1 support 
this prediction. The interfacial tensiois, ycg contains three tertns, 
which, in principle, can all be controlled. We performed exper- 
iments with fixed cell-cell adhesion and vatying gel properties. 
For the purposes of this study, we used iV-cadherin transfected 
CHO cells. Theit adlsesive properties were quantitatively as- 
sessed by measuring aggiegate surface tension (for details on the 
method, see refs. 12, 24, 33, and 34). We sised gels with differing 
chemical composition. The relative importance of cell-cell and 
ceii-rostrix interactions has been investigated quantitatively (11). 

Results in Fig. 1 / and J show that collagen at concentration 
of 1.7 mg/mi .is analogous to a permissive scaft'old with small 
ycg/ET. Collagen at concentration of 1.0 and 1.2 mg/ml atid the 
RGD containing neurogel match more the definition of the 
nor.perniissive gel with Isigii ycg/Ej. Ti.ese gels favor much less 
(collagen), or not at all (neurogel), the dispersion of the ceils into 
the scaffold, tiuis facilitating fusion. During our measurements, 
we did not observe the collapse of the fused rings that, according 
to the model predictions, should eventually take place (Fig. 2). 

The influence of Scaffold on 3D Structure Formation. The above 
results demonstrate the well known fact that scaffold properties 
affect cellular structure (6, 11). The specific mechanism of how 

the gel influences pattern evolution depends on its detailed 
chemistry, and in general is not easy to discern. 

Cells exert traction forces on their substrates and surroutjding 
3D matrices (35-38). In the case of collagen, numerous sliidies 
of this phenomenon have been pet fortned (39-43). Tlte evolu- 
tion process in this case has some striiiing features. As the fusion 
of aggregates takes place, the pattern noticeably contracts, at 
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Fig, 3. Time evolution of the ceiiular pattern InFig. If. Notetlie initial strong 
contraction (see also Fig. 4). Also note that there are many more ceils migrat- 
ing into tl-ie gei tiian thh im.^ge rriiglit suggest. Becasjse the fluorescent signal 

of individual ceiis is much weai<;er than that of the fused aggregates, they are 
hard to spot (however, see Fig. 5). The collagen layers above and under the 
ring are too thin to observe any appreciable cell migration out of the plane of 
the ring. 



least until ^60 h (Figs. 3 and 4). For higher collagen concen- 
tration, contraction is more drarnatic, whereas la the case oi the 
neurogel no similar effects were observed (results not shown). 
Ccmiraction re.suits from the C1I.O ceils puiting on collagen 
fibers. Fuithermore, as aggregates fuse, at some point the 
pattern assumes a starburst appearance (see spikes at ^-60 h, but 
not before, in Fig. .5 for the 1.0 mg/ml collagen gel). Cells extend 
outward from the ring suggesting that by this time a network of 
radially aligned collagen fibers has developed in the vicinity of 
the aggregates. (Foi' an explicit visualizatiot! of such an ariay, see 
ref. 42.) It is remarkable that no similar effect is observed inside 
the ring even at 144 h (Fig. .5). The probable reason collagen 
fibers do not align inside the ring is that the vectorial sum of the 
isoJropicaily acting liaciion forces is zero. (An analogous can- 
cellation makes the electric field in the interior of a conducting 
spherical shell to be zero.) 

In Fig. 4, we quantified contraction in teriiis of the total area 
defined by the outer perimeter of the ring. The curves are 
exponential fits to the data in i lic toi iny4exp(-i/rcg) -t- B {A and 
B are constants). The quantity Teg defines a characteristic time 
scale of contraction. (For the 1.0 mg/mi collagen gel, the one we 
will use in later analysis. Teg 57 h.) 

Kiiietlcs of Aggregate Fusion. The transpatency of collagen allows 
optically following structure evolution. Fig. 5 shows the time 
variation of the boundary between two adjacent aggregates in 
the 1.0 mg/ml collagen gel. A measure of fusion is the instan- 
taneotis value of the angle formed by the two aggregates. As 
aggregates coalesce, the angle between the tangents to their 
boundaries (drawn from the point where they join) approaches 
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Fig. 4. Cel!-matrl.x interaction induced contraction in collagen gels. The 
figure shows the change of the entire are.^i inside the outer perimeter of the 
rings in Fig. 3. Squares, circles, and diamonds stand for collagen concentration 
1 .0, T.2, and 1 .7 mg/ml, respectively. Curves are exponential fits to the data. 
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Fig. 5. Time course of aggregate fusion for collagen concentration 1 .0 
mg/ml. Bright field images of the cellular pattern were acquired with a X4 
objective. {Top) Evolution of the cellular boundary between adjacent aggre- 
gates. Note the radial, exclusively outward (top) directed migration of ceils at 
around 60 h. (Middle) The entire pattern at selected times {/.eft, 35 h; Right, 
144 h). (Bottom) variation of the angle between aggregate boundaries in Top 
as function of time. The curve is an exponential fit to the data (see text for 
details). 



180^ The curve in Fig. S is an exponential fit to the daia in ilio 
form C[\ -exp(-/,/Tcc)l (C is consiaot), witii Tqc, 23 ii. Here, 
the quantity Tcc defines a time scale of aggregate fusion. 

Dsscisssion 

We have manually printed a simple, but nontiivial structure, a 
ring of spherical aggregates, each containing thousands of ceils 

with specific adhesive propeities. We have shown, both exper- 
imentaliy and in cotnpuie! simulations, that under appropriate 
conditions the initially coniigtious aggregates, positioned along 
a circle, fuse int o a toroidal construct. If aggregates were printed 
in multiple layers, tliey would presumably fuse in both the 
horizontal and vertical directions, thus forming a iumenous 
organ-like module. 

The major outcome of this wotk is the det)ionstratio!i that 
spherical cell aggregates can be used as building blocks in tissue 
engineering applications. In particular, with the development of 
autotsialed, co!n[)uter-aided disjsensers ot bioprinters (which 
already exist for the delivery of cells; refs. 13 and 14) they could 
be used in the capacity of bio-ink. 

The biophysical basis for bio-ink is tissue liquidity, a concept 
proposed by Steinberg (18). The differential adltesion hypothesis 
(DAH) provides the molecular foundation for tissue liquidity 
(18, 19). DAfl explains the liquid-like behavior of ceil popula- 
tions in terjns of tissue surface and interfacial tensions, gener- 
ated by adhesive and cohesive interactions between the compo- 
nent subunits. In vitro experiments have demoiisliated that tissue 
surface tension is a well defined physical parameter, wliiclt 
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characterizes ihe eqiiilibriuni shape of miiiticeiiiilar aggregates. 
Measured values of the tensions in many cases account for the 
observed rautuai envelopment behavior of tissues (12, 34). In 
vivo experiments give further support to DAIf (44, 45). 

According to FJAH, the lowest energy configuration of an}; 
tissue fragmeni (con tain ing motile ceiis, utiifoi'mly acihesive tive i 
their entire surfaces) embedded in a medium to which it adheres 
weakly, is a sphere, Tiius, the final pattern in the 1.0 and 1.2 
mg/nii collagen gels and the neurogel should also be a single 
spherical aggregate. However, as the simulations demonstrate, 
the system may be trapped in long-lived metastable intermediate 
states (Fig. 2), which correspond to particular liised conforma- 
tions. This may provide sufficient time to dissolve the gel, thus 
freezing the desirable configuration and transferring the result- 
ing tissue construct into a bioieactor for maintetiatice. 

The cell aggregates we used in this work contain cells of only one 
ty[>e. Coiriplex organs contain several ceil types. Tliere ase indica- 
tions that using aggregates comprised of heteroceUuIar populations 
wouid lead to nontriviai structures. It has been shown, boiJi 
experimentally (12, 34) and in computer simulations (26, 27), that 
paitern evoluiion in sijrting is consistent with the predictions of 
DAH, and is indeed controlled and driven by interfacial tension. 
Moreover, when the mixture is composed of ceils erf tissues that 
are neighbors in normal development, in the course of sorting 
they recover their physiologica! cotifigiu'ation (46, 48). 

Liquidity has been emphasized for embryonic tissues because 
it is primarily at this stage of development that tissues and cells 
must actively move to eventually give rise to organs. Here we 
used aggregates of A'-cadherin transfected CHO ceUs, whose 
liquidity has beets established earlier (49). It i.s quite likely that 
when structure formation is induced by the methods de.scribed 
here, iiggicgatcs of otiici cell types (in particular .stem cells) will 
behave similarly to embryonic tissu'.;s, as far as their liijuidity and 
ability to fuse is concerned. 

Spherical cell aggregates can be made only of adhering cells. 
Fibroblasts, for example, normally do not adhere to each other 
directly, only tbivH.gh ihu extracellular matrix. Such cells can 
either be temporarily genetically manipulated to expjess cell 
adhesion molecules (50, 51) or embedded in a population of 
adhering cells. 

Delivering or printing ceil aggiegaics (instead of individual 
cells) ititc) 3D scaffolds offers several advantages. Aggregates are 
prebuilt small tissue blocks, thus their fusion immediately results 
in 3D structiires. Because lliey contain many thousands of cells, 
printing titne could be dramatically reduced and cell survival 
greatly improved. Properly designed composite aggregates (with 
more than one cell type and potentially cotstaits ing extracellulat 
matrix; i.e., "multicolor bio-ink") can enhance the creation oi 
desired complex tissue constructs. Finally, the inescaf/able haish 
mechanical conditions encountered in presently available print- 
ers are more critical for cells thais cell aggregates. [.Tet-based ceil 
printers (15), due to their small orifice, are presently not 
appropriate for dispen.sing aggregates. Devices with the poten- 
tial to print aggregates of several hundred microns (14) use 
micropipettes as cartridges with presstire-operated extrude:rs, 
and thus provide more gentle conditions.] 

Our work is based on the analysis of cell- hydroge! interac- 
tions in 3D. The computer simulations suggest that for the 
ftision of the aggregates to take place strict conditions on the 
embedding gel must be imposed. Our experiments support this 
prediction. In the case of collagen, because of the tendency of 
cells to reorganize the matrix, the situation is complex, and the 
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2. Langei, R. & Vacanti, .1. P. (1993) .Vcfewcc 260, 920-926. 

3. Bonassar, L. .1. & Vacami, C. A. (1998) / Cell Biochem- Suppl 30-31, 
297-303, 

4. L'j'saght, M. J., Ngujf, N. A. & Sullivan, K. (1998) Tissue Eng. 4, 231-238. 



interpretation of the experiments in terms of only gel- 
aggregate iiiterfacial tensioii igiiores important details. Figs. 4 
and 5 imply that pattern evolution involves two distinct time 
scales. Teg 57 h and Tcc ^ 23 h for the 1.0 tng/ml coliageii 
matrix, which characterize respectively the dynatuical aspects 
of cell-gel interaction (i.e., how long it takes for the cells to 
reorganize the matrix) and cell-cell interactions (i.e., how long 
it takes for the aggregates to fuse). It is the competition 
between these interjictions that determines the final pattern. 
Because for 1.0 mg/ml collagen Tcg/Vcc '■"■='2.5, contraction is 
considerably slower than fusion, and the ring has time to 
stabilize. For 1.7 mg/ml collagen concentration, corttraction 
dominates pattern evolution and ihe fused ring does not: have 
time lo stabilize. 

The ideal hydrogel for cell aggregate printing must allow 
cells to survive and provide favorable conditions for postprint- 
ing self-assembly. There already exist several candidates for 
such gels: thermo-reversible gels (52, 53), photo-sensitive gels 
(54-56), pH-sensitive gels (57), and gels sensitive to specific 
molecular entities (58, 59), Detailed studies must be per- 
formed lo identify the optimal candidate for .specific cells or 
their mixtures. 

Our mode; to describe the observed patter;) cvoiiitiori clearly 
represents a sttcing (iveisitnplificatio;) of the coHtplex nature of 
both the cellular and scaifoid system, lis major limitation is that 
it assumes tJiat cells are identical point particles. Thus, in 
particular, it ignores the change in their shape, a factor known 
to be important in cell-matrk interact hm. Simiiat iy, :t heats the 
embedding gel as a system without internal sti lictutc, which is 
not an aciequate characterization of a network of fibers, in 
particular the otte sf^attnirig a collaget) ttiattix. Furtherntore, 
cellular rearrangement within the model is exclusively driven by 
energy minimization (as implemented l)y the Monte Carlo 
method), thus, no information can be gained on the true time 
evolution of tlie pat tern. The entire process is conlxoUed by a 
single parameter, the aggregate-gel interfacial tension. It is 
intriguing that, despite these limitations, the similarity betw'een 
the sturralated and experimentally obtained shapes is so strong. 
This .suggests that our model may incorporate the most saHent 
features of the observed phenomenon. Possible improvement, of 
the model includes the treatment of cells as extended objects, 
capable of changing their shape. This can be accomplished for 
example within the large state Potts model (26, 27). Following 
the true time evolution of the pattern would require the detailed 
characterization of forces between cells and gel, and necessitate 
molecular dynamics based simulations instead of the Monte 
Carlo approach. 

In conclusion, we have demonstrated that closely placed cell 
aggregates in 3D gels can self-organize into metastable tis.sue 
constructs of desired shape. It was also shown that adhesive and 
mechaniccil properties of the embedding hydrogels are critical 
for cell aggregate fusion. We have developed a mathensatical 
model that accurately describes our experimental findings. 
When exteiided to incorporate more tealistic cell and hydrogel 
properties, it may be used for the design and optimization of 
scaffold properties. Finally, on I he basis of our findiiigs, we have 
proposed that cell aggregates could be used as self-assembling 
bio-ink in automated delivery or printing devices. 
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ABSTRACT 

Self-assembly is a fundamental process, which drives structural organization in both 
inanimate and living systems. It is in the course of self-assembly of cells and tissues in 
early development that the organism and its parts eventually acquire their final shape. 
Even though developmental patterning through self-assembly is under strict genetic 
control it is clear that physical mechanisms must underlie the formation of complex 
structures. Here we show, both experimentally and using computer simulations how 
tissue liquidity can be employed to build tissue constructs of prescribed geometry in 
vitro. Spherical aggregates, which consist of many thousand cells and form due to tissue 
liquidity were implanted contiguously into biocompatible hydrogels in circular geometry. 
Depending on the properties of the gel, upon incubation the aggregates either fused into 
toroidal three-dimensional structure or their constituent cells dispersed into the 
surrounding matrix. The model simulations, which reproduced the experimentally 
observed shapes indicate that the control parameter of structure evolution is the 
aggregate-gel interfacial tension. The model-based analysis also revealed that the 
observed toroidal structure represents a metastable state of the cellular system whose 
lifetime depends on the magnitude of cell -cell and cell-matrix interactions. Thus these 
constructs can be made long-lived. We suggest that spherical aggregates composed of 
organ-specific cells may be used as "bioink" in the evolving technology of organ 
printing. 
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Self-assembly is the fundamental process, which generates structural organization across 
scales (1). Histogenesis and organogenesis are examples of self-assembly processes in 
which through cell-cell and cell-extra cellular matrix interactions the developing 
organism and its parts gradually acquire their final shape. In the present work we use both 
experimental and computational approaches to demonstrate how the self-organizing 
properties of cells and tissues, the basis for morphogenesis, can be exploited to build 
three-dimensional (3D) biological structures of prescribed geometry. 

Tissue engineering (2-7) offers the opportunity to study self-assembling processes 
during histo- and organogenesis in vitro under controlled conditions. Tissue engineering 
aims to not only creating desirable organs, but also to better understand the fundamental 
mechanisms and principles of biological organization in general and morphogenesis in 
particular. Classical tissue engineering is based on seeding cells into biodegradable 
polymer scaffolds or gels, culturing and expanding them in bioreactors for several weeks 
and finally implanting the resulting tissue into the patient's organism where the 
maturation of the new organ takes place. 

It has recently been suggested to use cell aggregates, instead of individual cells as 
building blocks in tissue engineering (8,9). Cell aggregates have traditionally been used 
as a powerful tool to understand the principles of cell-cell (10) and cell-matrix adhesion 
(11) as well as cell sorting (12). In addition, rapid prototyping technology has 
successfully been applied for computer-aided deposition of cells in gels in order to create 
3D tissue constructs (13,14). We introduce the concept of "bioink", spherical cell 
aggregates, which upon implantation or "printing" into a scaffold ("biopaper") have the 
ability to fuse into 3D organ structures (1 5-17). 
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The ability of cell aggregates to fuse is based on the notion of tissue fluidity 
(18,19), according to which embryonic tissues in many respect can be considered as 
liquids. In particular, in suspension or on non-adhesive surfaces, various multicellular 
aggregates round up into spherical shape similarly to liquid droplets (12). We 
hypothesize that closely placed aggregates in appropriately chosen 3D gels can fuse to 
form tissue constructs of desired geometry. 

To demonstrate the feasibility of such a proposition we used aggregates of 
genetically transformed cells with controlled adhesive properties and arranged them to 
form a ring in gels of different chemical and mechanical properties. Our results 
demonstrate that contiguous aggregates under appropriate conditions, defined by the 
composition of the embedding gel indeed can fuse into structures of specified 
morphology. We have also constructed and experimentally validated a mathematical 
model of cell aggregate fusion with strong predictive power. 

Taken together these results suggest that cell aggregates can successfully be used 
as building blocks in tissue engineering technologies. They also provide support for the 
concept of self-assembling bioink and thus justification for the use of cell aggregates in 
the evolving organ printing tools. 

Materials and Methods 

Cell Aggregate Preparation. Chinese Hamster Ovary (CHO) cells transfected with 
N-cadherin (courtesy of A. Bershadsky, Weizmann Institute, Rehovot, Israel) were 
infected with histon binding H2B-YFP retrovirus (kindly provided by R.D. Lansford, 
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Beckman Institute at California Institute of Technology). Confluent cell cultures (3-4x10^ 
cells/75 cm^ TC dish) grown in Dulbecco's Modified Eagle Medium (DMEM, Gibco 
BRL Grand Island, NY; with 4500 mg/L glucose, L-Glutamine and pyridoxine 
hydrochloride and supplemented with 10% FBS (US Biotechnologies, Parkerford, PA), 
10|Lig/ml of penicillin, streptomycin, gentamicin, kanamycin, 400|ig/ml geneticin and 
MEM sodium pyruvate) were washed twice with Hanks Balanced Salt Solution (HBSS) 
containing 2mM CaCh, then treated for 10 minutes with Trypsin 0.1% (diluted from 
2.5% stock, Gibco BRL, Grand Island, NY). Depleted cells were centrifuged at 2500 
RPM for 4 minutes. The resulting pellet was transferred into capillary micropipettes of 
500 ^m diameter, incubated at 37 C" with 5% CO2 for 10 minutes. The firm cylinders of 
cells removed from the pipettes were cut into fragment (500 |xm height), then incubated 
in 10-ml tissue culture flasks (Bellco Glass, Vineland, NJ) with 3 ml DMEM on a 
gyratory shaker at 120 RPM with 5% CO2 at 37°C for 24-36 hours. This procedure 
reproducibly provides spherical aggregates of similar (-500 |j,m diameter) size. 

Cell Aggregate-Gel Structures. NeuroGel™ (a biocompatible porous poly[N- 
(hydroxypropil)methacrylamide] hydrogel) disks of 10mm diameter and 2 mm thickness, 
containing RGB fragments (kindly provided by Stephane Woerly, Organogel Canada, 
Quebec) were washed three times with DMEM to eliminate the storage medium. This gel 
has been shown to provide favorable conditions for spinal cord repair (20,21). A 0.5 mm 
wide, 0.5 mm deep circular groove was cut into a disk, then filled with 10 aggregates, 
placed contiguously. The groove was refilled with the gel to completely embed the 
aggregates. This structure was incubated at 37°C, 5% CO2 for 72 hours in a tissue culture 
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dish containing 10 ml DMEM, washed with PBS and finally embedded in Tissue-Tek® 
OCT Compound (Electron Microscopy Sciences, Fort Washington, PA). The structure 
was slowly cooled (1 C^/min) to -20 C" in a Nalgene freezing container (Nalgene 
Labware, Rochester, NY). To visualize aggregate fusion, at the end of the experiment 
cryosectioning was performed with a Reichert 2800N Frigocut cryotome (Reichert Jung, 
Amsberg, Germany), yielding 10-16 [im thin slices mounted on microscope slides. Slices 
were examined on an Olympus IX-70 inverted microscope with fluorescent attachment at 
4x magnification. 

To tune the strength of cell-gel interaction further fusion experiments were 
conducted in rat-tail collagen type I (Sigma-Aldrich, St. Louis, MO). Collagen was 
dissolved in IN acetic acid then treated with Ham's F12 medium with sodium 
bicarbonate. Under room temperature this mixture gels in a few minutes depending on 
concentration. The gel-aggregate structure was achieved by creating a ring of ten 
aggregates on the top of a previously (almost) solidified collagen layer, then covering 
with liquid collagen that embeds the aggregates after gelation. These samples were 
incubated under the same conditions as described above. Working with 1.0 mg/ml and 
1 .7 mg/ml collagen, the samples were transparent, thus it was possible to follow pattern 
evolution in time by phase contrast and fluorescent microscopy. 

Cell survivability was checked with Trypan Blue (Invitrogen, Carlsbad, CA) at 
the end of the fusion experiments. A minimal number of uniformly distributed dead cells 
were found. 
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Modeling structure formation via cell aggregate fusion. To investigate shape 
changes of the aggregate structure in gel, we constructed a three-dimensional model, in 
which the sites of a cubic lattice are considered to be occupied either by cells or by gel 
volume elements. The total interaction energy or the Hamiltonian of the system is written 
as 

H= E/(a„a,0, (1) 

<r, r'> 

where rand r' label lattice sites, and <r, r' > signifies summation over neighboring 
sites, each pair counted once. First, second and third nearest neighbors are included, and 
we assume that a cell interacts with the same strength with all the 26 cells it comes into 
contact with. To specify occupancy, we assign a spin value, g , to each lattice site with 
values 0 for a "gel particle" and 1 for a cell. The interaction energy of two neighbors, 
J(<j^,a^,), may take either of the values J{0,0) = -s^, J(l,l) = -s^^ or 

J(0, 1) = /(1,0) = -s^g . Here the positive parameters , and s^^ account for contact 

interaction strengths for cell-cell, gel-gel and cell-gel pairs, respectively. More 
specifically, these are mechanical works needed to disrupt the corresponding bonds. 

(Note that s^^ and s^^ are works of cohesion, whereas is work of adhesion per bond 

(22)). The strength of cell-cell interaction may be determined experimentally either 
directly (23) or by measuring the tissue surface tension (12,24), whereas cell-gel 
interaction is tunable via the concentration of RGD groups in the gel (6,25) or by the 
concentration of collagen. The gel-gel "bond energy" is an effective measure of gel 



7 



filament density, interaction and stiffness and is determined by the specific chemistry of 

the gel. 

The Hamiltonian in Eq. 1 may be rewritten by separating interfacial and bulk 
terms in the sum. As a result we obtain 

H = y^B^ + const. (2) 

Here B^^ is the total number of cell-gel bonds, and y^g = (s^^ +Sgg)/2-s^g is proportional 

to the cell-gel interfacial tension (22,26,27) or spreading coefficient (28,29). The 
remaining terms of the Hamiltonian do not change as the cellular pattern evolves. 

Our model is inspired by earlier efforts aiming at computer simulations of cell 
sorting (the morphogenetic phenomenon in which one of two, initially randomly 
intermixed cell populations sorts out and becomes surrounded by the other). These were 
mainly concerned with finding the appropriate motility rules (30). Sorting was 
successfully simulated using the large-Q Potts model (26,27). In our simulations, due to 
the need to handle large numbers of cells, shape changes of individual cells are presently 
not accounted for. However, as far as the evolution of the cellular pattern is concerned, 
we have adopted the same stochastic procedure used in the Potts model -based studies 
(26,27). 

The evolution of the system is followed using Monte Carlo simulations (31), 
relying on a random number generator of L'Ecuyer with Baym -Durham shuffle (32). The 
program identifies the cells on the aggregate-gel interface, picks one of them randomly, 
and exchanges it with an adjacent gel particle chosen by chance. The corresponding 
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change in adhesive energy, ISH , is calculated and the new configuration accepted with a 
probability P = \ if A//<0 or P = eKp{-j3 AH) if A//>0, where J3 = l/E^, is the 
inverse of the average biological fluctuation energy, analogous to the thermal fluctuation 
energy (33), k^T (Ar^-Boltzmann's constant T- absolute temperature). Whereas in 
statistical mechanics this energy characterizes thermal agitation in a system of atoms or 
molecules, in the case of cells it may be thought of as a measure of the spontaneous, 
cytoskeleton driven motion of cells, able to break adhesive bonds between neighbors via 
membrane ruffling (34). By definition, a Monte Carlo step (MCS) is completed when 
each cell in contact with the gel has been given the chance to move once. During each 
MCS the interfacial sites are selected in random order. The gel boundary is treated as a 
fixed physical limit of the system and cells are constrained to move within the gel. 

Results 

Structure formation: simulations. The result stated in Eq. 2 indicates that the 
evolution of the cellular pattern within the gel is governed by a single parameter, Jcg/Ej- , 

which, for cells with specific adhesive properties, is controlled by the chemistry of the 
gel. The theoretical analysis shows that once the characteristics of the gel are 
appropriately tuned, efficient fusion of adjacent aggregates takes place. A gel that 
promotes effective cell motility lowers the interfacial tension in comparison to the tissue 
surface tension (i.e. the interfacial tension between an aggregate and surrounding tissue 
culture medium). This is illustrated in the simulation shown in Fig. 1 A-B performed for a 
low interfacial tension (y^/Ej.= 0.25). The cells spread in the bulk of the gel and as a 

result of their high motility the shape of the cellular pattern evolves towards the lowest 
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energy state, which corresponds to a uniform distribution of cells within the gel. (The 
spheroidal structure in Fig. IB eventually disperses.) Under optimal cell-gel interface 
properties, expressed in our model by a certain range of the spreading coefficient, fusion 
occurs resulting in a three-dimensional toroidal structure. An example of this kind is 
depicted in Fig. II-F (y^/^^ = 0-9). 

Energetics during structure formation. Some shapes correspond to local minima 
of the interaction energy: they represent metastable configurations. These are identified 
from plateaus in the plot of the total interaction energy vs. MCS, and are important for 
organ engineering, for they can be made long-lived. The initial configuration in Fig. II 
progresses towards such a long-lived metastable state. This is illustrated in more detail in 
Fig. 2, where the energy of the metastable state is essentially unchanged in the entire 
interval between 10"* and 6 x 10"^ MCS. Eventually the fused toroidal state becomes 
unstable: a particular section of the structure becomes gradually thinner. This corresponds 
to a slight negative slope of the energy plot, as shown in Fig 2B, and at about 10^ MCS 
the torus brakes down. Subsequent massive rearrangements lead to a pronounced energy 
decrease while the system evolves into three rounded aggregates. These remain stable 
because large spatial separations hinder their fusion into a single large spheroid. Similar 
simulations showed that metastability depends both on system size and interaction 
strengths. 

Once the cellular structure reaches the metastable state, it can be stabilized by 
dissolving the supporting gel. In the simulations this corresponds to increasing the value 
of the aggregate-gel interfacial tension. Indeed, if in the simulations shown in Fig. 2A, 
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this quantity is changed to y^^/£'j.=2 anywhere in the plateau region, the energy 
remains constant as long as the simulation is run (results not shown). 

Experimental realization of three-dimensional structure. To study the 
feasibility of engineering three-dimensional tissue constructs of prescribed geometry, we 
have "manually printed" (i.e. embedded) aggregates of living cells into biocompatible 
gels following insight form computer simulations. As indicated by our model, the ability 
of aggregates to fuse depends on the mutual properties of the cells and the gel, as 
expressed by the parameter Ycg/^r • The results of our experiments, summarized in Fig. 1 

and 3, support this prediction. The spreading coefficient contains three terms, which, in 
principle, can all be controlled. We performed experiments with fixed cell -cell adhesion 
and varying gel properties. For the purposes of this study, we used N-cadherin transfected 
CHO cells. Adhesion between particular cells can quantitatively be assessed, by 
measuring the surface tension of tissues composed of these cells, as it was done in a 
number of cases (12, 24, 35). Gel properties were varied by using materials with differing 
chemical composition. The importance of cell -cell versus cell-matrix interactions has also 
been investigated quantitatively (1 1). 

Results in Fig. IC-D show that collagen at concentration of 1.7 mg/ml is 
analogous to a permissive scaffold with small y^jEj.. On the contrary, collagen at 
concentration of 1.0 mg/ml and the RGD containing neurogel match more the definition 
of the non-permissive gel with high y^J Ej, These gels do not favor the dispersion of the 
cells into the scaffold thus facilitating fusion, as can be seen in Figs. lE-F, I-J and Fig. 3. 
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During our measurements we did not observe the collapse of the fused rings that 
according to the model predictions should eventually happen (see Fig. 2). 

Mechanism of 3D structure formation in collagen. It is well known that cells 
exert traction forces on their surrounding matrices (36-39). This may affect structure 
evolution, and indeed as Fig. 3 shows the circular configuration markedly contracts at 
least up to a certain point. The mechanism of structure formation in collagen, as shown in 
Figs. 1 and 3, can be understood using recent quantitative results on cell -collagen 
interactions (40-43). In case of a single cell aggregate embedded in collagen cells 
reorganize the originally homogeneous and isotropic array of collagen fibers. They 
initially tow the matrix symmetrically from all directions and simultaneously start 
spreading along them. This results in cell treadmilling with little net cell migration. With 
time, collagen movement slows down, whereas cell spreading away from the aggregate 
accelerates. When two aggregates are embedded into collagen, as a consequence of 
traction collagen fibers align parallel to each other and form a strap between the 
aggregates with subsequent cell migration along the strap (43). 

The structure evolution shown in Fig. 3, although more complicated, is consistent 
with the earlier findings on cell-matrix interactions. In the configuration of 10 aggregates 
(Fig. 3), each pair along the diagonal of the ring would give rise to a collagen strap in the 
absence of the others. Due to symmetry, the cumulative effect of all the aggregates is that 
no stripes form inside the ring (collagen fibers "do not know" between which aggregates 
to align). At the same time each aggregate can exert traction on the fibers outside of the 
ring. As fibers are towed they push on the aggregates. The pattern effectively contracts. 
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since no counterbalancing structure inside the circle exists. Using the time evolution of 
the pattern as showed in Fig. 3 (and a similar sequence for the higher collagen 
concentration, not shown) we measured the area defined by the fluorescent ring as 
function of time. Results are shown in Fig. 4. As the movement of collagen fibers outside 
the ring slows down for collagen at concentration of 1.0 mg/ml, the rate of contraction 
decreases rapidly, but remains almost steady for the higher collagen concentration. The 
less concentrated collagen preparation represents a thin gel on which the cells' ability to 
pull and thus migrate is limited (38). If we assume that area contraction is predominantly 

due to collagen fiber movement, the rate of contraction ^[7ir^(^)]= C(0 determined 

from Fig. 4 allows to obtain a measure for v(t), the instantaneous speed of collagen fiber 
movement by v(t) = C{t)/2'Kr{t). This speed decreases from its initial value of 0.02 
|j/min and 0.04 [i/min to 0.003 |i/min and 0.02 |j/min, respectively for collagen 
concentrations 1.0 mg/ml and 1.7 mg/ml. These results are consistent with the findings of 
Sawhney and Howard (43) (these authors used fibroblasts in their studies). 

Kinetics of aggregate fusion. Since the transparency of collagen allows to follow 
structure evolution in time, we were able to quantify the development of the fused ring as 
shown in Fig. 3. Fig. 5 A shows the time evolution of the boundary between two adjacent 
aggregates in the 1.0 mg/ml collagen gel. Note that up to 60 hours, there is minimal 
migration of cells into collagen and when migration is pronounced it occurs only 
outward, but not inward. This may be the indication that no alignment of collagen fibers 
inside the ring has taken place. A measure of fusion is the angle formed by the two 
aggregates, whose time dependence is shown in Fig. 5B. As the perimeter of the ring 
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flattens by the coalescence of aggregates, the angle between their boundaries approaches 
180 degrees. The steady increase of this angle shows that cells, at least at the aggregate- 
gel interface were alive. (Cell viability for the entire construct was checked as explained 
in the Materials and Methods.) 

Discussion 

We have manually printed a simple but non-trivial structure, a ring of spherical 
aggregates, each containing thousands of living cells with specific adhesive properties. 
We have shown, both in computer simulations and experimentally, that the initially 
contiguous aggregates positioned along a circle, under appropriate conditions fuse into a 
toroidal tissue construct. If aggregates were printed in multiple layers (layer-by-layer 
deposition) they would presumably fuse in both the horizontal and vertical directions, 
thus forming a vessel-like structure. 

The major outcome of this work is the demonstration that spherical cell 
aggregates can be considered as bioink in the sense that they have the capacity i) to be 
delivered by computer-aided jet-based or automatic cell dispensor-based "printing" and 
ii) to consolidate to form "text" or self-assembled histological or tissue constructs. 

The biophysical basis for bioink is provided by tissue liquidity. The liquidity of 
embryonic tissues was proposed by Steinberg (18) as a possible mechanism underlying 
early developmental patterning. The Differential Adhesion Hypothesis (DAH; 18) 
provides the molecular foundation for tissue liquidity. DAH explains the liquid-like 
behavior of cell populations in terms of tissue surface and interfacial tensions, generated 
by adhesive and cohesive interactions between the component subunits. In vitro 
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experiments have demonstrated that tissue surface tension is a well-defined intensive 
physical parameter, which characterizes the equilibrium shape of multicellular aggregates 
and the measured values of these tensions account for the observed mutual envelopment 
behavior of the corresponding tissues (12). In vivo experiments give further support to 
DAH (44-45). 

According to DAH the lowest energy configuration of any tissue fragment 
containing motile cells embedded in a medium to which it adheres weakly, consistent 
with the notion of tissue liquidity is a sphere. Thus, the final pattern in the 1.0 mg/ml 
collagen gel and the neurogel should also be a single spherical aggregate. However, as 
the simulations demonstrate, the system may be trapped in long-lived metastable 
intermediate states (Fig. 2), which correspond to particular fused conformations. This 
may provide sufficient time to dissolve the gel thus freezing the desirable metastable 
configuration and transferring the resulting tissue construct into the bioreactor for 
maintenance. 

Our model and cellular aggregates are clearly oversimplified representations of 
reality. The model assumes that cells are identical particles, which do not change shape in 
the course of pattern evolution, which is driven by energy minimization and is controlled 
by a single parameter, the aggregate-gel interfacial tension. The strong analogy between 
the stimulated and experimentally obtained shapes nevertheless strongly suggests that our 
model incorporates the most salient features of the observed structure evolution. 

The cell aggregates we used in this work contain cells of only one type; they 
allow printing only with a single "color". Complex organs contain several cell types. 
There are many indications that the use of multicolor bioink would lead to highly 
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nontrivial organ structures. It has been shown by many investigators, both experimentally 
(12,34,46-48) and using computer simulations (26,27) that sorting follows the behavior of 
immiscible fluids and is indeed controlled and driven by the interfacial tension. 
Moreover, when the mixture is composed of cells of tissues that are neighbors in normal 
development, in the course of sorting they recover their physiological configuration 
(46,47). 

Tissue liquidity has been emphasized for embryonic tissues since it is primarily in 
this stage of development that tissues and cells must actively move to eventually give rise 
to organs. Here we used aggregates of N-cadherin trasfected CHO cells, whose liquidity 
has been established (49). It is quite likely that when structure formation is induced by 
the methods described here aggregates of other cells types (in particular stem cells) will 
behave similarly to embryonic tissues, as far as their liquidity and ability to fuse is 
concerned. 

Spherical cell aggregates can be made of only adhering cells. Fibroblasts, for 
example, normally do not adhere to each other directly only through the ECM. Such cells 
can either be temporarily genetically manipulated to express cell adhesion molecules 
(50,51) or embedded in a population of adhering cells. 

Printing cell aggregates instead of individual cells offers several advantageous. 
First, aggregates are pre-built small tissue blocks, thus their fusion immediately results in 
3D structures. Second, since they contain many thousands of cells printing time could be 
dramatically reduced and thus cell survival greatly improved. Third, properly designed 
composite aggregates (more than one cell type and potentially containing ECM, i.e. 
multicolor bioink) can enhance the creation of desired complex tissue constructs. Finally, 
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the inescapable harsh mechanical conditions encountered in the nozzle of the printer are 
more critical for cells than cell aggregates. 

The real challenge for organ printing is not so much placing cells in 3D space, 
which has already been demonstrated to be technically feasible (13-17), but rather the 
capacity of these precisely placed cell aggregates to fuse and form histotypic and 
organotypic structures and then consolidate into mechanically stable and maintainable 
constructs. The latter is probably not possible without some sort of smart hydrogel. These 
gels must allow cell aggregates to fuse into the desired geometry. Our data strongly 
indicate that this can be accomplished only by assuring the right balance between cell- 
cell and cell-matrix (hydrogel) adhesion, as shown by Ryan et al. (1 1) for cell aggregate 
spreading on 2D substrates. 

Our work is based on the analysis of cell-hydrogel interactions in 3D. The 
computer simulations indicate that for the fusion of the aggregates to take place strict 
conditions on the embedding gel must be imposed. Our experiments show that by tuning 
the properties of the hydrogel, the fusion of aggregates made of N-cadherin transfected 
CHO cells could be altered. Collagen type I at concentration of 1.7 mg/ml does not favor 
aggregate fusion because it is too easy for the cells to migrate into the gel (the interfacial 
tension, in the simulations is too small). On the other hand collagen at concentration 

1.0 mg/ml or a special neurogel with optimal concentration of RGD fragments (6) allow 
Only restricted motion of cells away from the aggregates and thus promote fusion (y^^ is 
sufficiently large). 

The ideal hydrogel for organ printing must allow cells to survive and must 
provide favorable conditions for post-printing tissue fusion and self-assembly. There 
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already exist several candidates for such gels: thermo-reversible gels (52,53), photo- 
sensitive gels (54-56), pH- sensitive gels (57) and gels sensitive to specific molecular 
entities (58,59). Detailed studies must be performed to identify the optimal candidate for 
specific cells or their mixtures. The capacity to predict cell aggregate "fusogenic" 
behavior and guide the choice of the optimal gel are important attributes of the developed 
mathematical model. 

Several hypotheses have been proposed to explain the locomotive behavior of 
cells in 3D gels. According to one concept (60) it is the pore size of the gel that 
determines cell migration. However there is strong evidence that the presence in 3D gels 
of adhesive sites such as RGD fragments, especially in optimal concentration and 
arrangement is critical for cell attachment and motion (6, 61-64). 

Finally, the mechanical properties of gels are also critical for cell migration (36- 
43, 65,66). It would be beneficial for tissue engineering applications if the numerous gels 
now in use could be analyzed in terms of a few relevant physical parameters, as is the 
case (to an extreme degree) in our simplified model where all matrix properties (both 
chemical and mechanical) are reduced to a single quantity, the interfacial tension. It is 
clear, further specially designed experiments are necessary to elucidate this important 
point. 

In conclusion, we have demonstrated that closely placed cell aggregates in 3D 
gels can self-organize into metastable tissue constructs of desired form. It was also shown 
that adhesive and mechanical properties of the embedding hydrogels are critical for both 
cell aggregate fusion and the maintenance of the printed construct. We have developed a 
mathematical model that can be used for the design and optimization of hydrogel 
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properties and the prediction of the emerging biological structure. Finally, we have 
shown both experimentally and theoretically that cell aggregates could be used as self- 
assembling bioink for evolving organ printing technology. 

This work was supported by research grants from NSF (G.F.) and NASA (G.F., V.M., 
R.M.). 
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Figure Captions 

Fig. 1. Initial (upper row) and final (lower row) cell aggregate configurations in the 
simulations and in the various biocompatible gels. Panels A-B and E-F correspond to 
simulations with yagl^i - 0-25 and "Hcgj^T - 0-9, respectively. The ten aggregates, each 
containing 925 cells are one cell diameter from each other in the starting configurations. 
The finai configurations shown in B and F are reached after 25,000 and 50,000 MCS, 
respectively. Panels C-D, G-H and I-J correspond to CHO cell aggregates embedded in 
collagen gel of concentration 1.7 mg/ml and 1.0 mg/ml and a neurogel with RGD 
fragments, respectively. The nuclei of the cells are fluorescently labeled (see Materials 
and Mehtods) and the images of the cellular patterns were acquired with a 4x objective. 
The average diameter of the aggregates is 500 |i. 

Fig. 2. Metastable configuration at y^^ lE^ =1.1; (A) Fusion of 16 contiguous aggregates 

(123 cells each) takes place within the first 5000 MCS. The plateau of the energy vs. 
MCS corresponds to the thin ring conformation with an average height of 5 cell 
diameters. (B) The evolution of the energy during 5x10^ MCS and representative shapes. 
The final, rounded aggregates remained stable, and the energy practically constant up to 
10^ MCS (not shown). 

Fig. 3. Time evolution of the cellular pattern in Fig. IC. Note the initial strong 
contraction (see Fig. 4). Also note there are many more cells migrating into the gel than 
this image might suggest. Since the fluorescent signal of individual cells is much weaker 

than that of the fused aggregates they are hard to spot (see however Fig. 5). 
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Fig. 4. Cell-matrix interaction induced contraction in collagen gels. The figure shows the 
change of the entire area inside the outer perimeter of the rings in Fig. 3. Since few cells 
migrate into the gel, their effect on the area has been ignored. Squares and diamonds 
stand for collagen concentration 1.0 mg/ml and 1.7 mg/ml, respectively. 

Fig. 5. Time course of aggregate fiision for collagen concentration 1.0 mg/ml. Upper 
panel. The evolution of the cellular boundary between adjacent aggregates. Bright field 
images were acquired with a 4x objective. Lower panel. Variation of the angle between 
aggregate boundaries as function of time. 
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